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The nonlinear propagation of electron-acoustic solitary structures is investigated 
in a plasma containing k distributed (superthermal) electrons. Different types of 
localized structures are shown to exist. The occurrence of modulational instability 
is investigated. 



I. INTRODUCTION 



Superthermal particles in laboratory, space and astrophysical plasmas are often modeled 



v „,,,„„,,„„„ ^moQIt,,,.,,,,, „ T „ 



deviation from a Maxwellian distribution (the latter is recovered for infinite k). Our twofold 
aim here is to investigate the effect of superthermality on electrostatic solitary waves, and 
also on self-modulated wavepackets. 

Electron-acoustic waves (EAW) occur in plasmas containing two distinct temperature 



a a. 



These are 



electron populations (here referred to as "cold" and "hot" electrons) 
high frequency electrostatic electron oscillations, where the restoring force comes from the 
hot electrons pressure and the cold electrons provide the inertia jj Q], while ions plainly 
provide a neutralizing background. The phase speed v p h of the EAW is much larger than 
the thermal speeds of both cold electrons and ions but much smaller than the cold electrons 
(i.e., Vph }C , v p h,i <C v P h <C Vnhh). EAWs survive Landau damping in the region Th/T c > 10 
and 0.25 < n c0 /n h0 < 4 j^yj], where we have defined the temperature (T c , T h ) and density 
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FIG. 1: Variation of the pseudopotential ^f((j>) with <f> (left); the electric potential eft vs. £ (right). 
We have considered various values of k, and a = 0.01, j3 = 1, and M = 1. 



(n c , rift) of the electron constituents ('c' for cold, 'h' for hot). 



II. ELECTRON FLUID MODEL 



We consider a three component plasma consisting of inertial ("cool") electrons, k- 
distributed ("hot") electrons and stationary background ions. In a ID geometry, the dy- 
namics of the cold electrons is governed by the following normalized equations: 

On d(nu) du du dej) a dP 

dt dx dt dx dx n dx ' 

— +u— + 3P— = 0, — 1 = - (/3 + 1 )+n + /3 1-— V —\ 
at ox ox ox 1 \ k — | / 

We have scaled all relevant physical quantities as: cold electron density n = n c /n c0 ; fluid 

speed u = u c /v ; electric potential = $/$o! time t = tu pc ; space x = x/\ ; pressure 

P = P/n c0 k B T c ; we have defined: v = (k B T h /m e ) 1/2 , A = (k B T h /4:7in c0 e 2 ) 1/2 and oj~} = 

(47m c0 e 2 /m e )~ 1 / 2 , and also the density- and temperature- ratio(s): (3 = n ht0 /n C:0 and a = 

T c /T h . 



III. ARBITRARY AMPLITUDE SOLITARY EXCITATIONS 



Anticipating stationary profile localized excitations, we shift from variables {x, t} to x = 
x — Mt, where M is the solitary wave speed, scaled by v (defined above). We obtain 
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FIG. 2: Soliton existence region (Mi < M < M2) for different temperature ratio a values, versus 
f3 for k = 3 (left panel); versus n for j3 = 1 (right panel). 



u = M^l--^j,u = M- (M 2 +20 - 3n 2 a + 3a) 1/2 , and P = n 3 . Poisson's equation thus 



leads to a pseudo-energy balance equation 



where the "Sagdeev" pseudopotential function \l/(0) reads [6| 
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Soliton Existence 



In order for solitons to exist, we need to impose 



'(0 = 0) = and 



0) < 



(where the prime denotes differentiation wrt 0), leading to the (true sound speed) threshold 



Mi = 
ment 



K-3/2 



3a 
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i F 2 {M) 
value; \l/(0) is real for 

12 



1 / 2 

. An upper limit for M is obtained by imposing the reality require- 
^(•/OU^max > (where max is a limit on the electrostatic potential 
< 0; see Fig. [T^i). The region thus obtained is depicted in Figure 
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FIG. 3: Envelope type solitary excitations: bright type (left panel) and dark type (right panel). 

IV. MODULATED ELECTRON- ACOUSTIC WAVEPACKETS 

We consider small (e <C 1) deviations of all state variables, say S (= n,u,(f)), from the 
equilibrium state, viz. S = + 'E^L 1 e n Y^i=- n S^e il ^ kx ~ ui \ and allow for a weak space- 
/time-dependence of the l-th. harmonic amplitudes S^ nl \ In what follows, we ignore the 
pressure term (cold electron model) for simplicity, and set a = n^o/n^o (= 

The 1st order (~ e 1 ) expressions provide the EAW dispersion relation u 2 = fc a^_" ; along 
with the amplitudes of the first harmonics. The 2nd and Oth harmonics are obtained at 
order e 2 . Annihilation of secular terms at 3rd order yields a nonlinear Schrodinger (NLS) 
type equation: 

i£ + pg + «W* = 0. (3, 

where the amplitude ip = <jyp ((, r) depends on ( = e(x — v g t), r = e 2 t, while v g = % = 
and P and Q are dispersion and nonlinear coefficients respectively. 



A. Modulational instability 

Adopting standard procedure Q , we investigate the occurrence of modulational instabil- 
ity by considering a harmonic solution of (jSJ) and then a harmonic amplitude perturbation 
with (wavenumber, frequency) = (fc, cD). A nonlinear dispersion relation is thus obtained: 
Cj 2 = P 2 k 2 (k 2 - 2§|^i j0 | 2 ). Provided that PQ > 0, wavenumbers above k cr = (^r) 1/2 |^i,o| 
lead to (amplitude) modulational instability (MI). For PQ < 0, wavepackets are stable. 
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FIG. 4: PQ = (or k = k cr ) contours vs carrier wavenumber k and superthermality parameter k, 
or density ratio a. Left panel: a = 0.25 for the green curve; 1 for magenta; 2.5 for red; and 4 for 
blue. Right panel: k = 3 for green; 4 for magenta; 8 for red; and 100 for blue. 




FIG. 5: Modulational instability growth rate (normalized by its value for infinite k) versus pertur- 
bation wavenumber. Left panel: k = 100, 7, 5, 3.5 (top to bottom) for a = 0.5, k = 3.2. Right 
panel: a = 0.5, 1, 2, 4 (top to bottom) for k = 4.5, k = 7. 



B. Envelope solitons 

In the modulationally stable region (PQ < 0, in fact essentially for large wavelengths) 
S o litonS may OCCUr , i.e. exact Solutl on S U1 ft e fo rm & , = W , - ^[(f - 
Vr)/L]}. On the other hand, for > 0, "bright" envelope solitons occur in the form Q: 
^ = ?/'osech[(^ — Vt)/L\. In the above, is the asymptotic electric potential amplitude 
value, V is the propagation speed and L is the soliton width, while the positive constant d 
regulates the depth of the void (d — 1 for black solitons or d < 1 for grey ones). 
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V. SUMMARY 

Stronger superthermality leads to higher amplitude solitary excitations (as suggested by 
Fig. [1]). Both the cold electron temperature and concentration significantly effect on the 
soliton existence domain, as the upper Mach number limit M2 increases for higher "cold" 
electron temperature, while the sonic threshold Mi is decreased for higher n c o (see Fig. 
H|). The modulational instability growth rate may be reduced due to stronger superther- 
mality (see Fig. (5k), while (somewhat counter-intuitively) the presence of more excess hot 
(superthermal) electrons increases the instability growth rate (see Fig. [5b)- 
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